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ABSTRACT 

We use a combination of deep, high angular resolution imaging data from 
the CDFS (HST/ACS GOODS survey) and ground based near-IR Kg images to 
derive the evolution of the galaxy major merger rate in the redshift range 0.2 < 
z < 1.2. We select galaxies on the sole basis of their J-band rest-frame, absolute 
magnitude, which is a good tracer of the stellar mass. We find steep evolution 
with redshift, with the merger rate cx (1 + ^)3.43±o.49 optically selected pairs, 
and oc (1 + for pairs selected in the near-IR. Our result is unlikely to 

be affected by luminosity evolution which is relatively modest when using rest- 
frame J band selection. The apparently more rapid evolution that we find in the 
visible is likely caused by biases relating to incompleteness and spatial resolution 
affecting the ground based near IR photometry, underestimating pair counts at 
higher redshifts in the near-IR. The major merger rate was ~5.6 times higher at 
z ~ 1.2 than at the current epoch. Overall 41%x(0.5 Gyr/r) of all galaxies with 
Mj < —19.5 have undergone a major merger in the last ~ 8 Gyr, where r is the 
merger timescale. Interestingly, we find no effect on the derived major merger 
rate due to the presence of the large scale structure at z = 0.735 in the CDFS. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: interactions 
- galaxies: statistics 



1. Introduction 



Galaxy mergers are believed to be the chief mechanism driving galaxy evolution within 
the hierarchical framework. Although mergers are rare at the current epoch, the hierarchical 
framework predicts that the merger rate must have been higher at earlier epochs. Despite 
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its importance in understanding galaxy evolution, the quantification of the galaxy major 
merger rate an d its evolution with redshift is still an ill constrained and hotly debated issue. 
Patton et al. (119971 ) have derived a clear increase in the merger fraction with redshift till 
z ~ 0.33 with a power law index of 2.8 ± 0.9. However, extension of this work to higher 
redshifts has been riddled with controversy. Le Fevre et al. (120001 ) reported steep evolution 
of major merger rate till z ~ 1.0, with a power law index of 3.2 ± 0.6 using pair-counting in 
the optical band for identifying the merger candidates. Bundy et al. (12004 ) reported a much 
more modest evolutio n of m erger rate using K' band images for identifying major merger 
candidates. Lin et al. (120041 ) also r eporte d very weak evolution in the merger rates using the 
DEEP2 redshift survey. Bell et al. (120061 ) reported a fairly rapid evolution in merger fraction 
of massive galaxies between z ~ 0.8 and the cu rrent epoch, using the technique of projected 
2-point correlation function. Lotz et al. (120061 ) . on the other hand, reported that the major 
merger fraction remains roughly constant at ~ 7% ± 2% till z ~ 1.2, using a nonparametric 
technique for quantifying galaxy morphology. Some of the discrepancies in the results quoted 
above may be attributed to differences in sample selection criteria, different techniques to 
derive the merger rate and different definitions of major mergers used by various people. 

We have used a combination of high resolution imaging data from the HST/ACS 
GOODS survey and imaging data from the VLT foUowup of the GOODS-South field 
to quantify the major merger rate of galaxies and its redshift evolution using the technique 
of pair counting. Section [2] lists the datasets we have used in this work, and briefly explains 
the methodology that we have employed in this paper for deriving the merger rate. Section [3] 
details the sample selection criterion employed by us and the possible biases in our sample. 
Section H] explains the details of the photometry performed by us in the filter. Section [5] 
explains in detail how we identify major pairs of galaxies, along with correcting for possible 
contamination and incompleteness issues. Section [6] compares the differences in pairs iden- 
tified in the visible and the Near-IR. Section [7] deals with the calculation of major merger 
rate from our identified pairs and it's evolution with redshift. We conclude by discussing the 
implications of our results obtained in this paper in Section |8l We adopt a cosmology with 
ifo=70kmsec-iMpc-\ nM = 0.3 and nA = 0.7. 



2. The Data 

We have used version vl.O of the reduced, calibrated images of the Chandra Deep Field 
South (CDFS) acquired with HST/ACS as part of the GOODS survey (Giavalisco et al. 
2004). The SExtractor (Bertin and Arnouts 1996) based version rl.l of the ACS multi-band 
source catalogs was used for object identification. Spectroscopic redshifts were taken from 
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the redshift catalog of the VVDS (Le Fevre et ah 2004), G00DS/F0RS2 redshift survey 
(Vanzella et al. 2005, 2006) & from the IMAGES survey (Ravikumar et al. 2007). The near- 
m 3 k Ks band imaging data of the GOODS/CDFS region from the ESQ GOODS/EIS 
Release Version 1.5 (Vandame et al., in preparation) was used. 

2.1. The methodology 

The methodology that we employ in this paper is quite straightforward. First, we 
identify sources with known spectroscopic redshifts in the GOODS-S field. We then shortlist 
those galaxies with redshift in the range 0.2 < z < 1.2 and rest frame absolute magnitude 
brighter than Mj{AB) = —19.5 to obtain a volume limited sample of primary galaxies. We 
then identify neighbours within a projected radius of 5/i^qo kpc < r < 20/i]~oo kpc around 
each primary galaxy. If the difference in the apparent magnitude 6m of the primary galaxy 
and the neighbour satisfy the condition —1.5 < 6m < 1.5 in either the z or the Kg filter, 
the galaxy pair is designated to be a major pair in that particular filter. This is explained 
in Section [51 

We do not have the redshift for the secondary member of a major pair in most cases, 
and as such some of them can simply be foreground/background superpositions. Such a 
contamination is statistically corrected for by using number counts of objects in the field in 
the z or Kg filter and is explained in detail in Section 15.11 Our primary galaxy sample is 
then binned into three redshift bins. The fraction of galaxies undergoing major merger in 
any given redshift bin is then inferred from the fraction of galaxies existing in major pairs 
in the same bin, after applying a correction for foreground/background contamination. The 
evolution in the major merger fraction is then determined by looking for any change in the 
major merger fraction over the three redshift bins that we have. 

3. Sample selection and possible biases 

We cross-correlated the HST/ACS source catalog with the three redshift catalogs men- 
tioned above to yield objects with spectroscopic redshifts in the range 0.2 < z < 1.2. Since 
we have used spectroscopic redshifts from various surveys in the CDFS for the sample se- 
lection, our sample is liable to suffer from the same biases as those induced at the time 
of target selection by the spectroscopic surveys. The IMAGES survey is biased towards 
brighter galaxies with strong emission lines (Ravikumar et al. 2007). The shorter integra- 
tion times used in the IMAGES survey results in poor SNR for fainter galaxies and results 
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in difficulties in estimating redshifts for galaxies which exhibit only absorption lines or weak 
emission lines. Hence they are biased against faint, red early type galaxies. However, this 
bias is recovered to a large extent by the sample selection criterion of the G00DS/F0RS2 
redshift survey(Vanzella et al. 2005), which is tailor-made to exploit the red throughput and 
sensitivity of F0RS2. This preferentially selects red, faint early types galaxies, and provides 
a nice complimentary redshift c atalog to the IMAGES survey. Indeed this has been discussed 



at length by Ravikumar et al. (120071 ) (see their section 4.1) who have tested successfully the 
representativeness of the combination of IMAGES and F0RS2 sample. The VVDS survey, 
on the other hand, is quite unbiased owing to the simple selection criterion of apparent mag- 
nitude Iab < 24.0 and does not bias for or against any particular type of object (Le Fevre 
et al. 2004). The recovery success rate for the three redshifts surveys vary from 88% for the 
VVDS, 77% for G00DS/F0RS2 and 76% for the IMAGES survey. In addition to this, the 
above mentioned recovery rates are a function of magnitude of the object (see Le Fevre et 
al. 2004, Vanzella et al. 2005 & Ravikumar et al. 2007). Therefore, as in any flux limited 
survey, fainter objects with weaker or no emission lines (eg. red absorption line galaxies) are 
underrepresented. This caveat must be kept in mind while interpreting the results presented 
in this paper. Inspite of the limitations imposed by the flux-limited nature of the surveys, a 
combination of these three redshift catalogs yields a sample of galaxies which is reasonably 
representative of the field galaxy population out to a redshift of ~1.2. 

Galaxies with rest frame absolute magnitude brighter than Mj{AB) = —19.5 have been 
selected as the primary galaxies, leaving us with 695 galaxies. The luminosity cut off was 
selected to ensure that the sample does not suffer from any incompleteness (Malmquist bias) 
in the highest redshift bin. Rest-frame J band is significantly less affected by sporadic star 
formation than bluer bands, and is more representative of the stellar mass content. Fig. [1] 
illustrates that the resulting sa mple is rather well representative of the Schechter luminosity 



function from Pozzetti et al. (120031 ) at mean redshift Zmean ~ 0.5 & 1.0. The Schechter 



luminosity function is parametrized as 

ct>{M) = 0-9210 0*iO°-^("+i)(^^*-^^) X e-'o-<-^*-''' 



where the values of a, M* and </>* from Pozzetti et al. (120031 ) are given in Fig. [H 



Fig. [2] shows pictorially the sample selection criterion applied by us in terms of redshift 
and luminosity cutoffs. 
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4. Ks band photometry 

Near-infrared imaging observations of the CDFS have been carried out, as part of 
GOODS, in J, H, & Ks bands using VLT/ISAAC. In the present work, we have used the 
ESO GOODS/EIS Release Version 1.5 (Vandame et al., in preparation). This data release 
includes 24 fully reduced VLT/ISAAC fields in J & Ks bands, covering 159.1 and 159.7 
arcmin^ of the GOODS/CDFS region, respectively. 

We have carried out photometry with SExtractor on the J & Ks images, using DE- 
TECT_THRESH=1.0 ; DETECT_MINAREA=5 and DEBLEND_MINCONT=0.0005 in both 
the filters. Deblending is a major source of worry when studying close pairs of objects. How- 
ever, we are essentially immune to this problem as our selection criteria demands that the 
fainter member of a pair have at least 25% of the flux of the other member to qualify as a 
major pair. Hence even if SExtractor ends up splitting a source into two (or more) owing to 
a low value of DEBLEND_MINCONT, it will still not get counted as a legitimate major pair. 
The combination of DETECT_THRESH, DETECT_MINAREA and DEBLEND_MINCONT 
has been optimised after several iterations of experimenting with different values for these 
parameters and subsequent visual examination after overlaying the derived catalog onto the 
J & Ks band images. The differential number counts obtained using our SExtractor catalog 
are given in Fig. [3l 

5. Finding pairs of galaxies 

We have used a simple prescription for flnding pairs of galaxies by identifying neighbours 
within a projected radius of 20/ij"Qo ^pc for each of the 695 main galaxies. For each neighbour, 
we find the difference in the magnitude 5m of the main galaxy and the neighbour in both z 
as well as Kg filter separately. If the condition —1.5 < 5m < 1.5 is satisfied in at least one 
of the two filters, the pair is designated to be a tentative major pair in that particular filter. 
This condition ensures that the fainter member of the pair has at least 25% of the z or Kg 
band fiux of the brighter member. It is important to note here that it is entirely possible 
for a neighbour to be brighter than the primary galaxy. It may indeed happen that a fainter 
galaxy within a pair has a spectroscopic measurement, while the brighter galaxy has none, 
owing to observing constraints of the redshift survey. Not including brighter neighbours can 
lead to an underestimation of the derived merger fraction. We find that not including major 
pairs in which the neighbour is brighter than the primary galaxy, we will underestimate the 
number of major pairs by 34% in the Ks filter and 24% in the z filter. The importance of 
including brighter companions in the definition of major pairs has been missed by earlier 
workers in the field. 
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The simple selection criterion employed by us yielded a sample of 162 and 226 tentative 
major pair candidates in the Kg filter, and in the z filter, respectively. We do not have 
the redshift for the secondary member of a pair in most cases, therefore some of them can 
simply be foreground/background superpositions. Such a contamination can be statistically 
corrected for by using number counts of objects in the field and is explained in section 15.11 
However, in the cases where the secondary member does have a redshift, we have used it 
to rule out those apparent pairs from the list which are obviously at different redshift as 
compared to the host galaxy. This left us with 125 and 181 major pair candidates in the Kg 
filter, and in the z filter, respectively. Furthermore, a lower limit of 5/;,7ookpc is imposed 
on the separation between the major pair candidates and the main galaxy in order to avoid 
confusing clumpy extensions of the main galaxy itself as a neighbour. This left us with 
119 and 170 major pair candidates in the filter, and in the z filter, respectively. These 
major pair candidates, are then segregated into three redshift bins and are listed in Tabled! 
Fig. m shows a montage of z band HST/ACS images of some of the confirmed major pairs, 
identified in the z filter, where the secondary member has the same spectroscopic redshift 
as the host galaxy. Similarly, Fig. [5] shows Kg band images of some of the confirmed major 
pairs identified in the Kg filter. Note that there is a one-to-one correspondence between the 
first four objects in Fig. H] and Fig. [5l 

We have used the z as well as the Kg filter to define major mergers because the former 
filter allows to find companions very near to the primary galaxy (thanks to the exquisite 
spatial resolution of the ACS camera) and the latter is less affected by star formation and 
more representative of the stellar mass (e.g. Bundy et al. 2004). 



Since we do not have redshift for the secondary member of a pair in most cases, some 
of them can simply be foreground/background superpositions. Such a contamination can be 
statistically corrected for by using number counts of objects in the field in the z or Kg filter to 
estimate the probability of finding an object within projected 5/i7oo kpc < r < 20/1]^qo kpc of 
the main galaxy and within ±1.5 magnitude of the main galaxy. This statistically expected 
number can then be subtracted off the number of neighbours observed around a given host 
galaxy to yield the number of true neigbours for that particular host. Mathematically, this 
statistical correction is calculated as 



5.1. Estimation of the foreground/background contamination 




'm+1.5 



(2) 



- 7- 



where, N(m') is the differential number counts of objects in a given filter, m is the 
apparent magnitude of the host galaxy in the same filter, and ^2011'^^ kpc '"s'^roo '^p'^ 
projected radius in degrees corresponding to 20h^QQ kpc and S/i^qq kpc at the redshift of the 
main galaxy. In the Kg band, these number counts were established by using the SExtractor 
catalog as explained in section |H In the z band, we have estimated the number counts using 
the catalog published by the GOODS team. 

Equ. [2] gives us the number of sources that are expected to be present within ±1.5 
magnitude and within projected S/i^oo^PC ^ r < 20/;.7ookpc of the host galaxy by pure 
chance coincidence. It may be noted that this statistical correction is only calculated (and 
subtracted) in those cases where the neighbor does not have a spectroscopic redshift. If the 
neighbor has a spectroscopic redshift, we know for sure if it is a real pair or not, hence no 
need to apply the statistical correction. 

This number is calculated for every host galaxy in our sample which has a neighbour 
without a redshift. The final number is listed in Table [H segregated into three redshift bins, 
summed over all the host galaxies (with neighbours sans redshift) in a given redshift bin. 



5.2. Test of our contamination calculation 

We tested the accuracy of our statistically estimated foreground/background contam- 
ination derived in Section IHTT] in two ways. First, we compute the ratio of the statistically 
expected neighbours to the number of major neighbours actually found in projection as listed 
in Table [TJ This can be termed as a sort of conta minatio n rate. This contamination rate is 



then compared to earlier works of Le Fevre et al. (120001 ). The contamination rate that we 
find here (in the z band) is 42% ± 6.3%. This is in excellent agreement with the contamina- 
tion rate of 49% ± 13% (as derived from Table 3 of Le Fevre et al. 2000). This agreement 
gives us additional confidence in our estimate of the foreground/background contamination. 

In addition to comparison with earlier works, we also checked our contamination rate of 
42% ± 6.3% with the contamination rate derived using the subset of our major pairs where 
we have spectroscopic redshifts for both members of the pair. There are 65 major pairs in the 
z band where both members of the pair have spectroscopic redshifts. In order to establish 
how many of these 65 major pairs are indeed at the same redshift, we plotted a histogram 
of Sz/{1 + z) = {zhost — Zneighj/iX + ^host) for thcse 65 major pairs. This is shown in Fig. [61 
It turns out that there are a number of pairs with —0.006 < 5z/{l + z) < 0.006, which are 
presumably located at the same redshift, with the residual difference in redshift 5z/{l + z) 
being attributed to redshift measurement errors and peculier velocities of the members. 
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A 5z/ {1 + z) = 0.006 corresponds to delta velocity of 1800km/sec. We were forced to use 
this seemingly large value of Sz/{1 + z) due to the errorbars on the spectroscopic redshifts 
of the objects. The quoted errorbars in the redshift measurements are of the order ~0.002 
for VVDS and G00DS/F0RS2 and -0.006 for the IMAGES survey. In addition to this, 
these quoted errorbars are to be treated as lower limits due to the way in which they are 
calculated (see eg. Vanzella et al. 2005). Since we are using three different redshift surveys, 
we have taken to using the worst errorbars i.e. Sz/{1 + z) = 0.006. This implies that the 
delta velocity of 1800 km/sec is to be treated as an upper limit, and the true delta velocity 
is liable to be smaller. 

In Fig. [6] many objects have 6z/{l + z) much larger than 0.006 (so much so that most of 
them are out of the frame in Fig. [H]). The members of such pairs are treated as having discor- 
dant redshifts and are essentially foreground/background superpositions. From this plot, we 
can then derive the fraction of pairs that are actually foreground/background superpositions 
(i.e. the contamination rate). This contamination rate is found to be ~ 69% ± 13%. On 
first sight this might appear larger than the contamination rate of 42% ± 6.3% that we have 
obtained using our statistically estimated foreground/background contamination (although 
notice the large errorbars since we are dealing with small number statistics here). However, 
we checked the spectra of the pairs having discordant redshifts ourselves. It turns out that 
in many cases, the redshift determination can be quite poor indeed, specially in the low ex- 
posure and low resolution VVDS (as parametrised by their redshift flag). So in cases where 
the neighbor has a VVDS redshift flag 2 (75% confidence) or worse, there is a much greater 
risk of the pair being designated as having discordant redshifts. In addition to this, one also 
has to allow for the fact that some of the pairs using redshifts from the other two redshift 
catalogs can also be erroneously counted as discordant (for the same reason as above). Such 
pairs account for ~15/45 discordant redshift pairs in the z band. This can significantly swing 
the number of pairs with discordant redshifts. If we discard these 15 discordant cases and 
recalculate the contamination rate, we get 60% ± 14% as the contamination rate, which is 
in much better agreement with the contamination rate of 42% ± 6.3% that we have obtained 
using our statistically estimated foreground/background contamination. 

In the nutshell, our contamination calculation is robust and the contamination rate com- 
pares well with the contamination rate reported in literature, as well as with the subsample 
of our major pairs where we have spectroscopic redshifts for both members of the pair. Also, 
as a result of our eagerness to throw out prospective pairs as having discordant redshifts due 
to spurious redshift determination in some cases, we are going to end up with an estimate 
of merger fraction which is a bit too low and thus it should be treated as a lower limit. 
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5.3. Source blending and photometric completeness 

A major bias when studying close pairs of objects may be caused by the fact that sources 
are blended when their separation is comparable to the PSF. It particularly affects ground 
based Kg images because the Kg PSF (~0.5") is much larger than the HST optical PSF 
(~0.1"). This is illustrated in Fig. [TJ where we have plotted the angular separation between 
the main galaxy and neighbours identified in the z filter against the z magnitude of the 
neighbour, segregated into three redshift bins. A neighbour here is identified as any galaxy 
within 5/i]~ookpc < r < 20/ij~Qokpc of the primary galaxy, regardless of it's magnitude. The 
resultant list of 1085 neighbours includes major as well as minor neighbours, and constitutes 
the pool from which the major pairs are identified. In Fig. [TJ the open circles denote 
neighbours which have Kg magnitudes, whereas the green crosses denote neighbours which 
do not have Kg magnitude. As is evidenced in Fig. [71 there are a large number (~43%) 
of such neighbours which do not have Ks band magnitudes. This has the consequence of 
reducing the pool of neighbours from which the major pairs are identified in the Ks filter 
leading to an underestimation of the major pair fraction in the Ks filter. Most of the 
neighbours which do not have measurable Kg magnitude are rather faint, with the green 
crosses clustering towards the faint end. Although the 5/i]~oo kpc inner cutoff guards against 
blending in the Ks filter at lower redshifts, in the highest redshift bin, the angular size 
corresponding to S/i^oo^PC is small enough to allow legitimate neighbours to be identified 
in the z band within < 1.0" of the host galaxy, where blending in the Kg filter can be a 
serious concern. On further investigation, the bias is indeed found to be more severe in case 
of close neighbours, with ~54% of the neighbours within 1" of the main galaxy not having 
measured Kg magnitude. Furthermore, even out of the neighbours within 1" of the main 
galaxy which have a Kg magnitude ascribed to them, ~17% are actually blended with the 
host galaxy. Incompleteness of the Kg band catalog due to detection limit and blending is 
therefore the main reason for missing out Kg magnitudes for many of the neighbours. This 
implies that the major merger fraction that we derive using the Kg filter is a lower limit. In 
addition, this effect is more severe at higher redshifts, leading to an underestimation of the 
rate at which the merger fraction evolves with redshift. 

In the z band on the other hand, we have verified using the z band number counts in the 
GOODS-S field, that the photometric catalog is complete down to the faintest magnitude 
expected of a neighbour for the faintest primary galaxy (Mj = —19.5), even at the highest 
redshift (~1.2). 
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6. Visible versus near-IR identification of pairs 



As reported by Bundy et al. (120041 ). the selection of pairs in the visible bands may be 
affected by star formation, i.e. some low mass companions may be brightened in visible, 
leading to an overestimate of the actual fraction of major pairs. In order to check this 
assertion, we compare the z — Ks color of the primary galaxy to that of the neighbours 
which were identified in Section [51 This is shown in Fig. [8] (left). Not surprisingly, in the 
Ks filter we preferentially select neighbours which are redder than the primary galaxy, while 
a selection in the z filter selects bluer neighbours. In the cases where the z — Ks color 
of the neighbour is roughly similar to that of the primary galaxy, they are picked out as 
major neighbours in both the filters. It must be noted however that neighbours are equally 
distributed on both sides of the [z — Ks)host={z — Ks)ne.i nhhnv.r hne. In particular, we do not 



find evidence for the claim made by Bundy et al. (120041 ) that satellite galaxies tend to be 



bluer than the primary galaxies. In our work the s econdary galaxies can be brighter than the 



primary galaxies, whereas the Bundy et al. (120041 ) companions were always fainter than the 



main galaxies. In order to avoid any incompatibility in comparing our results with those of 
Bundy et al., we replotted Fig. [8] (right) using only those neighbours which are fainter than 
the host galaxy. We again find that the neighbours are equally distributed on both sides of 
the {z — Ks)host={z — Ks)neighbour huc. This shows that intrinsically the neighbour can be 
both redder as well as bluer compared to the host galaxy. 



7. The pair statistics and the merger rate evolution with redshift 

Table [1] lists the primary galaxy/neighbour statistics from our work, segregated in three 
redshift bins for each of the two filters, z and Kg. We have defined the major merger fraction 
(i.e. the fraction of galaxies likely to merge) as 



rf \ -^pairs ~l~ -^pairsi^COnf -N^proj. „ _ /q\ 

J[^) = ^ U.O [6) 

^^host 



where Nhost is the number of primary galaxies in each redshift bin, Npairs is the number 
of galaxies existing in major pairs within projected S/i^qq kpc < r < 20/i]~qq kpc of the primary 
galaxies in either of the two filters, Npairs{conf.) is the number of galaxies spectroscopically 
confirmed to be existing in major pairs and Nproj. is the number of major neighbours expected 
to be found within 5/i]~ookpc < r < 20/i]~ookpc of the primary galaxies by pure chance 
coincidence, estimated statistically using number counts of sources in either of the two filters. 
The factor 0.5 is the fraction of close pairs that are likely to merge as estimated by Patton et 
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al. (119971 ) for z=0 galaxies. This factor is possibly a function of redshift (Lavery et al. 2004), 
but the form of this dependence is not well established, and is one of the more uncertain 
factors in our study. While we have corrected our pair fraction for the foreground/background 
contamination, we do not have at our disposal measurements of relative velocity for the 
two galaxies found in close pairs in most cases. For some of the pairs where we do have 
spectroscopic redshifts for both the members, we have put a cutoff on delta velocity At> < 
1800 km/sec for identifying true pairs. As explained in Section [5^ this Av < 1800 km/sec 
is to be treated as an upper limit owing to the errorbars on the spectroscopic redshifts used 
and the true delta velocity is liable t o be much smaller. This justifies our use of the factor 



0.5 in Equ. [3] following Patton et al. (120021 ). who have studied z< 0.55 galaxies in pairs with 
dynamical measurements, and find that roughly half of those with Av < 500 km/s, are likely 
to be merging systems. 

The evolution in the merger fraction f(z) is plotted in Fig. [H In the left panel, the 
dashed line is the best fit curve of the form f{z) = /(O) x (1 + 2;)° to our z filter datapoints 
(filled circles) from this work, whereas the solid line is the best fit curve to our Kg filter data- 
points (filled triangles). The lower redshift points taken from literature use sample selection 
criterion which are inconsiste nt wit h our work and are shown only for comparison. The open 



stars are from Patton et al. (Il997l ). the open squares from Patton et al. (l2000l ). (120021 ) and 



the open circle is from the Millennium Gala xy Ca talogue (De Propris et al. 2005). All these 



works with the exception of Patton et al. ( 119971 ) choose their primary sample in the rest 



frame B, —21 < Mb < —18. The difference in the slopes obtained from the two filters z 
and Kg can be traced back to the photometric incompleteness and blending issues in the Kg 
filter as explained in detail in Section 15. 3[ 

The merger fraction shows strong evolution with redshift in both the z as well as the 
Kg band. However, our derived evolution for merger fraction suffers acutely from a lack of 
consistent datapoints at redshift~0.0. The GOODS dataset is based on a relatively small 
solid angle area (as expected for a deep survey), and is not ideal for deriving the merger 
fraction at redshifts lower that ~0.2, given the small sample size in this redshift range. 
On the other hand, the lower redshift estimates of merger fraction that we have available 
from the literature are inconsistent with our sample selection criterion and are hence not 
comparable. It is for this reason that we utilize the derived parameters from the fit that we 
have obtained using only our dataset (Fig. [9](a)), even though it comes at the cost of larger 
error in terms of the merger rate that we finally derive. 

To get an idea of how much our result is likely to change if we had a pivotal datapoint 
at redshift ~0.0, we have refitted the power law curves, this time including the lower redshift 
points from literature in obtaining the fit. This is shown in Fig. [9] right panel, where the 
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dashed line is the best fit curve to our z filter datapoints (filled circles) "plus all the five 
lower redshift points from literature. Similarly, the solid line is the best fit curve to our Kg 
filter datapoints (filled triangle) plus the five lower redshift points. The tacit assumption 
here is that the hypothetical redshift ~0.0 datapoint is some sort of an average of all the 
low redshift datapoints reported in literature. We notice that even in this case, the merger 
fraction shows strong evolution with redshift in both the z as well as the Kg band, with the 
rate of evolution in the Kg filter being slightly shallower on account of the incompleteness 
and blending issues described earlier. Since this fit was done solely with the purpose of 
giving an idea to the reader of the sensitivity of our derived results to the presence/absence 
of a redshift~0.0 datapoint, the results from this fit are not used anywhere in our work. 



We used a prescription similar to that of Lin et al. (120041 ) to convert the merger fraction 
(fraction of galaxies likely to merge) into a comoving merger rate (i.e. number of merger 

events/ M-pc" j Gyr) 



Nr^g{z)=^.hn{z)f{z)r~^ (4) 
where is th e comoving number density of galaxies (obtained by integrating the 



Pozzetti et al. (120031 ) luminosity function over the luminosity range of interest), /(z) is the 
merger fraction derived above, the factor 0.5 converts the number of merging galaxies into 
number of merger events and r is the merger timescale (assumed to be 0.5 Gyr). This yielded 
a merger rate of 2.08 x 10~^ Mpc"'^ Gyr^^ at the current epoch, which evolves by a factor 
~ 5.6 to 1.16 X 10~^ Mpc"'^ Gyr^^ at z = 1.2 using Kg data. For the z band data, these 
numbers are 1.16 x lO^^Mpc^^ Gyr^^ at the current epoch, which evolves by a factor ~ 15 
to 1.73 X 10-3 Mpc^^ Gyr"^ at ^ = 1.2. Integrating over the timescale from z = 1.2 to 
the current epoch, we find that 41% ± 3.5%(45% ± 13% for the z filter) of galaxies with 
Mj < —19.5 have undergone a major merger in the last ~ 8 Gyr. The merger timescale is 
generally estimated from simulations to range from 0.1 to 1 Gyr (see Hernquist and Mihos 
1995), which brings an additional uncertainty to our result. For example if it were assumed 
to be 0.35 Gyr(Carlberg et al. 2000), it would lead to 59% as being the fraction of galaxies 
having experienced a major merger in the last ~ 8 Gyr. 



7.1. The effect of the Large Scale Structure in the CDFS 

As is well reported in literature (Le Fevre et al. 2004, Ravikumar et al. 2007), there are 
several large scale structures in the CDFS which show up as spikes in the redshift distribution 
obtained from various redshift surveys. The most prominent large scale structure is at a 
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redshift of 0.735. It represents a factor ~10 overdensity in terms of galaxy number density. 
This has the potential of significantly altering the results that we have obtained for merger 
fraction in the redshift bin 0.5 < z < 0.75. The detailed effect of such large scale structures 
on the cosmological relevance of the GOODS south field have been considered by Ravikumar 
et al. J2007h . 

In order to check the robustness of our derived merger fraction, we recalculated the 
merger fraction in the redshift bin 0.5 < z < 0.75 by excluding the sources embedded in the 
large scale structure at z=0.735 (defined as objects within 5v < 1500 km s~^ of z=0.735). 
This resulted in the exclusion of 70 sources from the bin 0.5 < z < 0.75 out of which 11 were 
major pairs (in the Ks band). The merger fraction dropped from 5.48% to 5.34% which is 
well within the reported errorbars. The derived merger rate is therefore not affected by the 
presence of the large scale structure. 



8. Discussion 

We have derived the major merger rates of galaxies in the redshift range 0.2 < z < 1.2 
using pair counting in both optical and near-IR bands. This work provides a robust estimate 
of the major merger rate upto redshift ~1.2 using a representative sample of near-IR selected 
galaxies, which are well within the photometric completeness of the source catalogs. We find 
steep evolution with redshift, with the merger fraction oc (i-)-2)3.43±o.49 f^j, optically selected 
pairs, and oc (1 + z)'^-'^^^'^-'^^ for pairs selected in the near-IR. The difference in the slopes 
obtained from the two filters z and Kg can be traced back to the photometric incompleteness 
and blending issues in the Kg filter. We find that the major merger rate evolves by a factor 
~ 5.6 from 2.08 x 10"^ Mpc"^ Gyr~^ at the current epoch to 1.16 x lO'^ Mpc"^ Gyr"^ at 
z = 1.2. This implies that between 41% and 59% of all galaxies with Mj < —19.5 have 
undergone a major merger in the last ~ 8 Gyr, assuming that the average timescale for a 
pair to merge is in the range 0.5 to 0.35 Gyr. 

Our result is in agreement with the recently published results by the COSMOS team 
(Kartaltepe et al. 2007) who report a power law index of 3.1 ± 0.1, as well a s other works 



such as Le fevre et al. fl2000f ). Carlberg et al. fll994f ). Patton et al. fl2002f ) (see Table 2 
in Kartaltepe et al. 2007). Reports of a significantly lower merger rate (Lin et al. 2004, 
Carlberg et al. 2000), might be due to the fact that the investigators have used a luminosity 
evolution model for sample selection. Indeed such a correction is not significant in our work 
since we have selected our sample in the rest frame near-IR. 



The disagreement in our Kg derived merger fraction evolution with that reported by 
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Bundy et al. (120041 ) might be due to their relatively small statistics (our sample size is ~ 4 
times larger than theirs), and that they did not account for neighbours that are brighter 
than the primary galaxy. Our study implies that the star formation enhancement of the 
companion does not affect severely the calculation of the merger rate. 

We notice that at least some of the discrepancies between results derived in this paper 
and those of other workers in the field can be at tributed to different definitions of major 



mergers. The methodology of Lotz et al. (120061 ) probes a merger at a different stage in 



the merging process (ongoing merger), whereas our definition of major merger is sensitive 
to pairs of galaxies in an early stage of merger (upcoming mergers) so that two distinct 
members are easily distinguishable. Studies with 3D spectroscopy are certainly better suited 
to identify on-going or post-merger stages. Such studies (Flores et al. 2006, Yang et al. 2007) 
have found a frequent presence of complex kinematics at z=0.6 (26% of Mj<-20.3 galaxies), 
which if related to mergers, is 4 to 5 times higher than paired/on goin g mer gers found in 



studies like this one (see also Lotz et al. 2006). Indeed Hammer et al. (120071 ) noticed that 
all studies of merger rate agree to find 5%±1% for the merger fraction at z=0.6. If complex 
kinematics are related to mergers, it may imply a long duration (1-2 Gyr) for the remnant 
phase during which gaseous velocity fields are severely perturbed or chaotic (using a simple 
ergodic argument). In the near future we will combine 3D spectroscopy of paired galaxies 
to identify the different time-scales of the merging process which might shape the galaxies 
as we observe them today. 

Finally we notice that although Toomre type (Toomre & Toomre 1972) tidal features 
are seen in close pairs of galaxies in the lower redshift bins, such features are conspicuously 
absent in candidate pairs at higher redshifts. This is likely due to the (1 + z)'"^ surface 
brightness dimming, which quickly pushes the low surface brightness tidal features below 
the detection limits at higher redshifts. 



We thank the Centre Franco-Indien pour la Promotion de la Recherche Avancee (CE- 
FIPRA) for financial assistance under project number 2804-1. A.R. would like to thank 
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Table 1. Pair Statistics 



z 


Nhost 


t\tF850LP 
pairs 




AJ-F850LP 
proj. 


jyKs 
pairs 


N^Jrsiconf.) 


jS^Ks 
proj. 


f(z)(F850LP) 


f(z)(Ks) 


0.20-0.50 


61 


4 





0.70 


5 





0.81 


2.70 ± 1.67 


3.43 ± 1.87 


0.50-0.75 


294 


48 


7 


19.56 


39 


3 


9.76 


6.03 ± 1.30 


5.48 ± 1.14 


0.75-1.20 


340 


98 


13 


42.94 


62 


10 


18.80 


10.01 ± 1.64 


7.82 ± 1.31 


TOTAL 


695 


150 


20 


63.20 


106 


13 


29.37 







Note. — Nfiost is the number of primary galaxies in each redshift bin. Npairs is the number of galaxies existing in major pairs 
within projected S/i^qq kpc < r < 20/ifQo ^pc of the primary galaxies in each of the two filters. Npairs{conf.) is the number of 
galaxies spectroscopically confirmed to be existing in major pairs. Nproj. is the number of major neighbours expected to be ' 
found within kpc < r < 20/i]~qq kpc of the primary galaxies by piire chance coincidence, estimated statistically using"^ 
number counts of sources in each of the two filters. Merger fraction is estimated as {Npairs + Npairs{conf .) — Nproj.) /Nhost x 0.5. 
The errorbars are estimated using Poisson statistics. 
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Fig. 1. — A comparison between the luminosity distribution of our sample galaxies in the 
redshift range 0.2 < z < 0.65 (left panel) & 0.65 < z < 1.2 (right panel ) resp ectively and 
the distribution predicted by the luminosity function of Pozzetti et al. (l2003l ). The solid 
line is the luminosity distribution predicted by the mean values of a, M* and 0* given by 
Pozzetti et al. fl2003ll . whereas the two dashed lines on either side demarcate the 1 sigma 
region around the mean as parametrized by the errorbars on a, M* and (f)* mentioned above. 
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VVDS + F0RS2 + IMAGES redshift completeness 
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Fig. 2. — A pictorial representation of the sample selection criterion applied by us in terms 
of redshift and luminosity cutoffs. Notice the prominent large scale structure at 2; ~ 0.735 
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Fig. 3. — The Ks band differential number counts obtained using tlie SExtractor catalog 
that we have constructed. 
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Fig. 4. — A montage of z band HST/ACS images of some of the confirmed major pairs, 
identified in the z filter, where the secondary member has the same spectroscopic redshift as 
the host galaxy. Each image is 40 x 40/i]"oo ^P^ iii size. The unique object ID of the host 
galaxy is marked on the top left of each image, the redshift of the host galaxy is marked on 
the top right and the redshift of the neighbour is marked next to it. The horizontal bar on 
the bottom right shows a scale of 1 arcsec. 
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Fig. 5. — A montage of Kg band images of some of the confirmed major pairs, identified in 
the Kg filter, where the secondary member has the same spectroscopic redshift as the host 
galaxy. Each image is 40 x 40/i|fQQ kpc in size. The unique object ID of the host galaxy 
is marked on the top left of each image, the redshift of the host galaxy is marked on the 
top right and the redshift of the neighbour is marked next to it. The horizontal bar on the 
bottom right shows a scale of 1 arcsec. 
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Fig. 6. — Histogram of {zhost ~~ Zneigh) /(I + Zhost) for the 65 major pairs where both members 
have a spectroscopic redshift. Note that only 20 pairs have —0.006 < 5z/{l + z) < 0.006 as 
seen in the figure. Most other pairs have 5z/{l + z) much larger than 0.006 so that they are 
out of the frame. 
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Fig. 7. — The angular separation between the main galaxy and the neighbour identified 
in the z filter Vs the z magnitude of the neighbour, segregated into three redshift bins as 
indicated. The open circles are neighbours which have measured Kg magnitude, whereas 
the green crosses are neighbours which do not have magnitude. The horizontal lines 
correspond to 20/ij~Qo ^pc and S/i^qo ^pc at the mean redshift of each bin. 
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z — K host z — K host 



Fig. 8. — The left figure shows z — Ks color of the primary galaxy Vs that of the major 
neighbour (—1.5 < 6m < 1.5) as described in Section El The open circles and dots represent 
neighbours identified in the z filter and in the filter, respectively. Not surprisingly, in 
the Kg filter we preferentially select neighbours which are redder than the primary galaxy, 
while a selection in the z filter selects bluer neighbours. The figure on the right is similar to 
that on the left except that only those major neighbours that are fainter than the primary 
galaxies are plotted. 
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Fig. 9. — The merger fraction evolution with redshift in z filte r (filled circles) and Kg filter 
(filled triangle s) . T he operi stars are from Patton et al. (119971 ). the open squares are from 
Patton et al. (120001 ). (120021 ) and the open circle is from the Millennium Galaxy Catalogue 
(De Propris et al. 2005). In the left figure, the dashed line is the best fit curve of the form 
/(z) = /(O) X (1 + z)" to our z filter datapoints from this work, whereas the solid line is 
the best fit curve to our Kg filter datapoints from this work. The lower redshift points from 
literature are just shown on this figure for comparison. In the right figure, the lower redshift 
points from literature are included in obtaining the fit so that the dashed line is the best fit 
curve to our z filter datapoints from this work plus all the five lower redshift points from 
literature, and the solid hne is the best fit curve to our Kg filter datapoints plus the five 
lower redshift points. 



